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[Abstract] A technique for incorporating distributed aerodynamic normal forces and 
aeroelastic coupling effects into a stability analysis model of a launch vehicle is presented. 

The formulation augments the linear state-space launch vehicle plant dynamics that are 
compactly derived as a system of coupled linear differential equations representing small 
angular and translational perturbations of the rigid body, nozzle, and sloshing propellant 
coupled with normal vibration of a set of orthogonal modes. The interaction of generalized 
forces due to aeroelastic coupling and thrust can be expressed as a set of augmenting 
non-diagonal stiffness and damping matrices in modal coordinates with no penalty on 
system order. While the eigenvalues of the structural response in the presence of thrust and 
aeroelastic forcing can be predicted at a given flight condition independent of the remaining 
degrees of freedom, the coupled model provides confidence in closed-loop stability in the 
presence of rigid-body, slosh, and actuator dynamics. Simulation results are presented that 
characterize the coupled dynamic response of the Ares I launch vehicle and the impact of 
aeroelasticity on control system stability margins. 
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Translational perturbation, trajectory frame 

Angular perturbation, trajectory frame 
State of f h elastic eigemnode 
Dynamic pressure 
Drag force, A 
Span reference area 

Aerodynamic moment arm of x th station 
Axial drag force coefficient 

Mode shape and slope of t mode at gimbal 
Mode shape and slope of i th mode at X th station 


I. Introduction 

A well-known technique for assessing the short-period dynamic stability of an aircraft or missile in atmospheric 
flight involves the use of a set of linearized equations of motion about some trimmed or nominal operating 
condition. For a large, flexible launch vehicle, the problem is less trivial in that many additional degrees of freedom 
must be considered to effect proper stabilization of the actual plant dynamics. Due to their necessarily low structural 
mass, large boosters are very flexible; sensed angles and rates are easily corrupted by vibrations and significant 
dynamic coupling can occur among the structural modes, propellant slosh modes, and actuator dynamics. 
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Recent work has been undertaken 5 6 to revisit the assumptions of the launch vehicle perturbation models 1-4 that 
have been used in control design studies since the beginning of the ballistic missile era. A traditional technique is a 
careful elimination of the dynamic coupling paths that are deemed negligible or too difficult to handle analytically 
via transfer function methods. 1 The advent of rapid matrix-based numerical computation and a state-space 
implementation of the equations of motion allow the construction of a high-fidelity, high-order linear model, the 
response of which can be readily assessed for the purposes of stability analysis and control design. 6 Such a model is 
dynamically complete in that all of the terms evolved from the linear expression of the kinetic and potential energies 
are accounted for in the coupled dynamics; that is, the system should conserve energy when integrated numerically 
in the absence of external forces. Conveniently, complex elastic cross-coupling effects that were intractable to 
implement on analog computers can now be incorporated into the model via elegant linear algebra operations. For 
typical launch vehicles that by design operate in a linear regime, a linear stability analysis model can accurately 
capture the actual physics without resorting to a complex and computationally intensive multibody/multiphysics 
simulation approach. 

The present work extends the previous results by investigating the effects of additional generalized forces due to 
coupling of distributed aerodynamic loads and lateral bending of the vehicle structure and illustrates a technique for 
incorporating the relevant expressions into the previously derived state-space dynamics model. While similar 
solutions for aircraft structures are well known in the aeroelastic community, a straightforward technique for 
quantifying impacts to closed-loop control of launch vehicles in the presence of a thrust force and additional degrees 
of freedom does not appear to be readily available. Since the thrust vector control (TVC) loop bandwidth is limited 
and the nozzle is commanded with respect to the thrust structure, thrust forces (follower forces) rotate and translate 
in inertial space due to vehicle bending and also affect the bending dynamics. Although the rigid-body dynamics 
may be unstable in the absence of closed-loop attitude control, the dynamic stability of the open-loop system in the 
presence of aerodynamic and thrust forces can also be investigated in order to quantify the conservatism of other 
tools (e.g., nonlinear simulations) for which the implementation of aeroelastic models is excessively complex or 
computationally prohibitive. 


II. Local Angles of Attack 

As discussed in [6], the coupled plant dynamics are desired in an explicit state-space form suitable for 
numerical integration or stability analysis; the system 

x = E~ 1 Ax + E~ 1 Bu ^ 

A B 

can be obtained from the system matrices A, B via an inversion after writing the system dynamics in the general 
descriptor form as 
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using the block diagonal approach, with the state vector 

x=[cp (p 2 n, n 2 ••• n k h , n 2 ... n k 5 sl s s2 ... 5 sn 5 sl s s2 ... . C) 

Here, A r represents the rigid-body dynamics in the pitch (yaw) plane, A f includes the diagonal matrix of 
eigenvalues and structural damping coefficients used to construct the homogeneous flex dynamics, and A s contains 
the slosh dynamics. We will use the approach of augmenting the A r and A / dynamics with the components of 
rigid and elastic angle of attack so as to capture both pitch damping and structural interaction due to distributed 
aerodynamics, 2 and the relevant input matrices B r , B f will be modified to include wind gust effects. 

In most missile stability analysis models appearing in the literature, aerodynamics effects are accounted for via 
rigid body normal force and moment coefficients that are reduced to an equivalent center of pressure and normal 
force slope with respect to a trajectory-referenced angle of attack (point mass aerodynamics). 3 The components of 
the rigid-body angle of attack in the plane of motion are given by small-angle linearization of the rigid-body pitch 
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angle and lateral velocity in a quasi-inertial frame which is fixed with respect to the longitudinal axis of the vehicle. 
At the center of mass, the rigid-body angle of attack is described by 


z w 

a=cp -y+j> 


( 4 ) 


where <p is the trajectory-referenced rigid-body perturbation angle, z is the lateral velocity normal to the 
trajectory, and w is the rigid body gust velocity. 3 A small angle linearization of the inverse tangent has been used, 
with z ,w<g:V . These assumptions are reasonable for a typical booster vehicle on a gravity turn trajectory. 

If the angle of attack is computed at a station on the vehicle other than at the center of mass, the rigid 
translational velocity relative to the freestream and the angular displacement affect the angle of attack, and the gust 
velocity is a local velocity. Local angular displacements and lateral velocities consist of the sum of rigid body 
motion and flexibility effects. 

The local angular displacement in the pitch plane of the x dl discrete section of the vehicle with respect to the 
freestream velocity, with x an integer varying from 1 to / , is given by 


Z 1 x • ™x V" / ' 1 . \ 

“x=v----<p+y+Z [Vi X n i --v ix n i ] ■ 


( 5 ) 


The local translational velocity of the X th station is the sum of the translational velocity of the center of mass, 
the translational velocity due to the rigid angular rate, and the elastic velocity of the associated node. 2 These terms 
comprise those with coefficients of V ~ l , the latter two of which are sometimes neglected in stability studies. The 
components of the local angle of attack are illustrated in Figure 1. 



-yVixHi 


Figure 1: Trajectory perturbation for an elastically deformed vehicle 


Since the linear elasticity model describes motion at distinct gridpoints, the vehicle line load data is converted 
from a continuous distribution into a set of discrete rigid partial derivatives. For each * station from 1 to / , there 
is a corresponding load vector 
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that is derived from piecewise integration of the continuous distribution. Each local coefficient does not take into 
account dynamic propagation of flow along the vehicle nor the effects of elastic motion or deformation on this flow: 
it is purely a static-load, rigid-body assessment of the normal force distribution. This distribution is a first-order 
discrete approximation of the true dynamic loads under the assumption that deflections are small and do not 
significantly affect downstream flow. 7 A typical discrete line load distribution and the continuous data from which 
it is derived appear in Figure 2. 


Aerodynamic Line Load Distribution (filtered) 



Figure 2: Typical discrete line load distribution 


The equations of motion for the flexible vehicle 6 can be augmented with the complete angle of attack expression 
in (5). The three fundamental coupled equations (rigid body rotation, rigid body translation, and lateral bending) 
can be written as 


I(p=qSY d l, 
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Nozzle inertial coupling, slosh dynamics, and thrust forces are included in the model 6 but not explicitly written 
here since no direct coupling from aerodynamics occurs via those terms. The right-hand sides of (7), (8), and (9) 
represent only aerodynamic forces resulting from rigid-body rotation, rigid-body translation, elastic deformation, 
and wind gusts. 
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III. Rigid Body Aerodynamic Damping 


In [6] it was shown that the rigid-body dynamics could be written as 
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°^N - c, j 
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da V a 
_ & C n q S 
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( 10 ) 


assuming only a rigid-body force and moment partial derivative for a acting at a center of pressure location l a . 
The requisite rigid coupling matrices must be constructed from the vectors 

/*=[/, i 2 ... i, f, (id 

which is an / X 1 column vector of x station displacements and 
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... i? r 


( 12 ) 


also an / X 1 column vector. The rigid-body matrices are modified with the equivalent rigid-body partial 
derivatives computed from (6), (1 1), and (12), such that 
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qSk, 


V 


(13) 


where k, is the sum of the elements of d x _ and 


k=l T J x , k=[l 2 Jd x 


(14) 


Note the pitch damping terms appearing in (13). An inspection of the state vector (3) reveals damping effects 
on both the translational and angular accelerations due to angular velocity of the rigid body. The matrix (13) 
accounts for the first three terms on the right-hand side of (7) and (8). 


IV. Aeroelastic - Rigid Dynamic Coupling 

As it will be required shortly, define the binary vector operator 

a*b=[a l b 1 a 2 b 2 ... a„b„ } (15) 

where a,b are row or column vectors of the same dimension, and the result is a row or column vector, accordingly, 
having as its elements the products of the elements of a ,b . This is equivalent to the 1 .* ' operation appearing in 
the matrix programming language MATLAB® (The MathWorks, Inc.). 

Let V ix , Vix be the mode shape and slope, or the elements of the 7 th mode eigenvector at the X th station. 
Define the gridpoint influence matrices 
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which are each of dimension kxl over k modes and / stations. The choice of gridpoints to represent the local 
deformation of the vehicle obviously should well-represent overall motion of the load structure (and hence the outer 
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mold line) and not local elasticity. The use of centerline gridpoints in this manner is generally only accurate for low 
frequency global deformations, such as the first and second bending modes. 

Rigid-body coupling from flex occurs due to force follower terms in addition to distributed aerodynamics. The 
rigid coupling matrix due to elastic force following is already given by 1 




-R(i g v g i+^O -^(^^2+^2) - -R(i g v g k+Vgk) 
Ry gl Ry' g2 ... Rq>' gk 




(17) 


In order to form each modal coefficient of the generalized displacement in (7), the relation can be written as the 
inner product of a row vector of sums over all stations and the generalized displacement of k modes, that is. 
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with the vector generalized displacement i? = [ / 7i ^2 • • • *7* ] • Using the operator (15) and the matrix in ( 1 6), the 

row vector in (18) can be rewritten as 
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and similar results follow to accommodate the remaining terms in (7) and (8), where 
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( 20 ) 
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It follows that the flex-to-rigid coupling matrix can be compactly expressed as 

[O,, 


( 21 ) 


( 22 ) 




qs\*j, r 


(23) 


Finally, the flex to rigid coupling is the sum of the force follower effects and aeroelastic coupling, such that 


^r//2 


(24) 


V. Rigid Excitation of Flex and Intermodal Coupling 

Both rigid and flexible motion contribute to the modal forcing functions. The excitation of the 7 th of k modes 
by forces resulting from displacements and velocities of all k modes is an modal coupling effect of otherwise 
orthogonal modes. The generalized force on the / th bending mode due to rigid and flexible aerodynamic loads is 
represented by the right-hand side of (9). 

f A typographical error appears in this expression in reference [6]. 
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A. Rigid Excitation 

In the absence of distributed aerodynamics, the rigid-to-flex coupling matrix is the null matrix due to 
orthogonality; 

^f/r = [® 2 *xi] • (25) 

The first three terms on the right-hand side of (9) can be partitioned and written as coefficients of the state 
vector such that 



and it follows from our previous developments that (20), (21), and (22) can be used to compactly rewrite (26) as 




qSyJx V'[ l x* d x) ~^-Wad x 


B. Bending Influence Coefficients 


Finally, the remaining two terms on the right-hand side of (9) must be incorporated into the flex dynamics. 
Assuming normalized eigenvectors, the unforced flexibility equations have the form 

( 28 ) 


[K,+ K,+ K. D,+D.\ 


The diagonal matrices K h ,D h are the stiffness and damping matrices for the second-order flex dynamics; K f 
is an augmented stiffness matrix resulting from elastic force following, 6 and K a > D a are the augmented stiffness 
and damping formed from the outer products of the mode shape and slope along the vehicle, weighted by the line 
load vector. Via expansion and reordering of summations, write the first vector of double sums in (9) as 
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Inspection of the matrix on the right-hand side of (30) reveals that this expression can be written as the sum of 
I kxk matrices. 
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Each x of I kxk matrices <x v is formed from the dyadic product of the .v lh columns of ip .V x multiplied 
by the line load scalar corresponding to the x 01 station, that is. 
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and finally, (30) can be written 
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The right-hand side of (32) is the product of a /c X /c matrix and a column vector of elastic displacements. From 
this, it can be generalized that the remaining tenn in the flexibility equation can be represented using 
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(34) 


With the aeroelastic coupling matrices defined, it is noted that are linear m dynamic pressure, reference 

area, and velocity, and that the matrix components are a function only of the flex eigenvectors and the line load 
distribution. Furthermore, for homogeneous stability of the elastic response in the presence of aeroelastic ity and 
thrust but no other generalized forces, it is sufficient that A f is Hurwitz, and 

s Jt[A 7 ]<0. (35) 

The well-known criterion (35) predicts aeroelastic stability of free- free vibration as a function of dynamic 
pressure and velocity, independent of the remaining dynamics and attitude control system. It is important to note, 
however, that while condition (35) is necessary for homogeneous stability of the bending dynamics in the absence of 
non-aerodynamic forces, the result does not generalize when the remaining dynamics and/or closed-loop attitude 
control are incorporated. The present model allows for the evaluation of the closed-loop eigenvalues in the presence 
of the fully coupled system dynamics, including propellant slosh, actuator servo-dynamics, and the full autopilot. 
Finally, the open-loop system from gimbal input to gimbal command is a single-input, single-output (SISO) system 
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that yields a frequency response. Closed-loop stability and stability margins can be readily determined using the 
Nyquist criterion. 


VI. Gust Response 

The rigid and elastic wind response of the vehicle is driven by local lateral velocity at each station along the 
vehicle centerline. The rigid response to a gust input is determined from the sum of the forces and moments due to 
gust inputs appearing at each x of / stations. The existing model inputs // = [/? $] are the nozzle position and 

acceleration from an actuator dynamics model. For gust penetration studies, the input can be expanded to include 
the unique lateral wind gust velocity at each station, resulting in an input vector 

u = [p P | Wj w 2 ... w,\ (36) 


and input distribution matrices B r> B f that are expanded to dimension 3x(2 + /) and 2kx(2 +/) , respectively. 
For typical ascent flight conditions, however, the forward velocity is sufficiently high that it is reasonable to assume 
that the vehicle encounters the change in wind velocity at every station simultaneously, hi this case, all w t =w as 
the gust velocity collapses to a single input. 

If a single input is assumed, it follows from the relationships derived above that the sum of the coefficients of 
over all * stations on the right-hand side of (7), (8) can be represented using the same rigid lateral velocity 
entries in (13). The rigid-body input coupling matrix is expanded such that 
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Likewise, the flex mput distribution matrix that accounts for the right-hand side of (9) is formed from the 
product of the mode shape matrix (16) and the load distribution (6), using the same coefficients as those for lateral 
velocity shown on the rightmost column of (26). The augmented flex input distribution matrix therefore becomes 
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( 38 ) 
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VII. Simulation 

The present formulation has been incorporated into a general missile autopilot stability analysis and design tool, 
FRACTAL. 6 FRACTAL allows the trajectory to be simulated at a particular flight condition under the assumptions 
of fixed (LTI) or time-varying (LTV) parameters so as to assess the impact of variations in a given vehicle 
configuration and autopilot design on the relative stability of the closed-loop system. The perturbation dynamics are 
inherently linear and the use of classical stability margins is appropriate. 

As is often the case with large flexible boosters, phase stabilization of the first bending mode is a requirement to 
ensure high robustness to bending mode uncertainty; as such, the system is conditionally stable in phase. Since 
large boosters tend to be aerodynamic ally unstable unless augmented with external aerodynamic stabilization 
devices, the system is conditionally stable in gain. Various assessments of coupled transient response and stability 
have been performed, a subset of which is discussed below. 
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A. Root loci of the aeroelastic dynamics 


Use of the flexibility equations (28) alone provides the open-loop eigenvalues of the elastic response in the 
presence of a thrust follower force and aeroelastic forces at various thrust levels and dynamic pressures and can be 
readily obtained by varying the parameters in (29). A representative configuration of the Ares I vehicle was 
incorporated into the present model and various root loci were computed at dynamic pressures ranging from zero to 
the maximum dispersed dynamic pressure for a given flight condition. 8 An example of the root locus of elastic 
eigenvalues is shown in Figure 3. 



Figure 3: Aeroelastic root locus near maximum dynamic pressure 


It is clear from the figure that the general trend is an increase in aerodynamic damping from the assumed 
structural damping value of Kbi~ 0.005 to as much as £*,-=0.01 for the lower frequency modes. Very little 
frequency shift (due to rotational coupling and force follower effects) is observed for this vehicle configuration; the 
dominant effect is due to the translational degrees of freedom. 

Conclusions regarding modal damping must consider the nature of the planar, rigid load approximation and the 
limited fidelity in the higher frequency modes. An increase in modal damping in the linear analysis is reassuring, 
but is not definitive evidence that modal damping is increasing. It would be unwise, therefore, to rely on this modal 
damping for control design. 


B. Stability margin assessment 

An investigation of stability margins using the integrated linear formulation allows one to determine the relative 
impact of aeroelasticity on the vehicle stability margins with all other dynamics (propellant slosh, nozzle degrees of 
freedom) included in the analysis. An LTI Monte Carlo analysis is conducted at each operating point of interest and 
various stability margins are extracted and verified using an automated process. The most significant impact to the 
open-loop stability with the autopilot included can be predicted from the preceding eigenvalue analysis: attenuation 
of flexibility due to increased structural damping has a marked impact on the elastic response. 

A comparison of the nominal frequency responses is presented in a Nichols plot in Figure 1 . The inclusion of 
the distributed line loads has the anticipated effect shown from the eigenvalue analysis. The amplitude of the first 
bending mode is decreased due to the increase in aerodynamic damping, resulting in an increase in high-frequency 
phase margin for the phase-stabilized first bending mode. Interestingly, slosh mode response is also affected due to 
coupling through bending. 
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phase 


Figure 4: Nichols response of coupled model 

For Monte Carlo evaluation of stability margins, several tens of thousands of cases are analyzed with various 
parameter dispersions sourced from differing vehicle trajectories. An example of a Monte Carlo analysis performed 
at various flight conditions is shown below in Figures 5 and 6. The effects of rigid line loads on the coupled 
open-loop autopilot stability are evident as a general trend of decreased aerodynamic gain margin near maximum 
dynamic pressure. Of course, this effect is well within the vehicle and autopilot design capabilities 5 , but it illustrates 
the more aggressive rigid body low-frequency aerodynamic response encountered via the use of rigid distributed 
loads data in the fully coupled model. 


Minimum Aero GM, 2000 Run Sets 


Minimum Flex AM 2000 Run Sets 




Figure 5: Monte Carlo minimum aerodynamic gain margin Figure 6: Monte Carlo peak elastic repsonse 


§ While Hie aerodynamic gain margin tends to decrease with the inclusion of rigid line loads, the rigid body phase margin 
increases due to aerodynamic pitch damping, potentially allowing for improved perfoimance via additional autopilot tuning. 
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The effects of aerodynamic damping on the higher frequency modes are presented in Figure 6. As was seen in 
the nominal stability analysis, the flex attenuation or minimum flex amplitude margin (typically driven by the 
second lateral bending mode) increases due to the aerodynamic damping effect. The trend throughout flight is an 
overall decrease in the open-loop elastic response sensed by the rate gyros of the attitude control system. Such a 
decrease is beneficial as decreased bending filter attenuation requirements potentially allow for unproved controller 
bandwidth and attitude tracking performance, given sufficient confidence in the aerodynamics model. 

C. Coupled transient response analysis 

A visualization tool has been constructed with a dynamics engine based upon the FRACTAL model 
architecture. The visualization system is intended to allow for rapid evaluation of autopilot response to gust loads 
and disturbances at various frozen-time flight conditions through animation. For example, the onset of a 
steady-state wind or wind profile can be simulated to evaluate the ability of the in-flight load relief system to turn 
the vehicle into the wind. Lateral bending, slosh, and nozzle deflections are depicted in real time, and trajectory- 
relative position and velocity errors are computed. An example of the output of the visualization tool is shown 
below in Figure 7. The response to a steady-state wind disturbance has been simulated to evaluate the bending 
dynamics and trajectory error. 


FRACTAL gust response 
Scale Factor 20x /40x /20x (angle / flex / slosh ) 



VIII. Discussion 

The present work has illustrated a technique for the incorporation of a planar aeroelastic model into the 
equations of motion for stability analysis of a traditional launch vehicle. The approach is based on the classical 
perturbation methods for stability assessment. While the model is linear and is a first-order approximation to the 
true dynamics encountered in a non-steady- state load condition, the general trends can be expected to hold so long 
as the vehicle is indeed fairly rotationally symmetric and the assumption of uncoupled motion in a single plane is 
realistic. While the present model cannot predict the complex interacting physics that can be simulated in coupled 
computational fluid dynamics (CFD)/fmite element model (FEM) codes, it benefits from being compact and 
extremely efficient to implement in a matrix-based programming language. Rapid execution allows for the 
formulation to be used within optimization routines and for rapid Monte Carlo stability analysis. Further work will 
compare the present results with evolving higher-fidelity data to assess the accuracy of the technique when 
compared with large-scale nonlinear multiphysics simulations of the vehicle dynamics. 
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From Science to Solutions 


Overview and Objectives 


• Motivations for a coupled aeroelastic model 

• Challenges in launch vehicle control 

• Existing approaches in analysis of flexible rocket dynamics 

• Assumptions and limitations 

• Model derivation and implementation 

• Example results 
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Motivations 


• From a vehicle integration perspective, aeroelasticity, loads, and closed-loop 
autopilot stability are considered decoupled problems 

Aeroelastic stability is verified via extensive analysis assuming free-free structural 
modes; loads are verified using resultant modal damping from active stabilization 
(controller in the loop) and aeroelastic damping (positive or negative) 

Control stability is verified using rigid-body aerodynamics with a coupled rigid-flex- 
slosh-actuator model 

• Recent work has focused on revisiting and improving heritage launch vehicle 
dynamics models using modern linear algebra tools 

Numerical recipes for constructing high order linear models have been developed 

These models are fully coupled within the linearity assumptions (well justified for 
traditional launch vehicle configurations) 

• Present work incorporates rigid line loads into the fully-coupled model 

Allow evaluation of autopilot stability in the presence of an aeroelastic system 
Enhance confidence of aeroelastic stability in the presence of closed loop control 
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Challenges in Launch Vehicle Control 


• Complex dynamic models 

• Nonlinear, coupled, flexible multi-body system with variable mass 

• Flexible body models can be very high order 

• Limited controller architecture 

• Linear controllers prefaced with bending filters are preferred based on flight heritage 

• Highly flexible 

• Sensed angles and rates are corrupted by vibration 

Structural oscillations must be mitigated to minimize bending loads, aerodynamics, and 
propellant motion 

• Conditionally stable (aerodynamics and phase-stable bending mode) 

• Large uncertainties in key parameters 

• Significant dynamic coupling (flex, slosh, aerodynamics, nozzle dynamics) 
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Perturbation Dynamics History 


• Small perturbations have classically been used in assessing short-period 
stability of aircraft about a trim condition 

Launch vehicle models were derived from similar assumptions and matured 
through the 1950s-1960s 

• Perturbation models are well suited to launch vehicles due to the small 
deflections from the trajectory encountered in the ascent environment 

• For example, gravity turn trajectories with small angles of attack 

Classical models were of reduced order, and derived in transfer function or block 
form for implementation on an analog computer or for computing frequency 
response 

Control design often performed using root-locus methods, resulting in a PD-type 
controller prefaced by bending filters that removed elastic response from the 
autopilot attitude and rate channels 
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Modeling Assumptions 


• General assumptions 




Vehicle is tetragonally symmetric and 
exhibits negligible pitch/yaw coupling 


Reference trajectory is a gravity turn where g 
can be ignored: 

i, _£siiiX c 
*c — 

Vehicle undergoes small perturbations from 
the reference trajectory, and elasticity can be 
represented in the linear range 



Propellant slosh can be modeled via a 
spring-mass-damper mechanical analogue 

Aerodynamics are well-behaved and can be 
linearized about zero angle of attack 



Quasi-inertial frame is a trajectory frame that 
translates axially with the vehicle (2-DoF) 
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Modeling Assumptions (II) 


• Vehicle is represented by a quasi-rigid base body on which orthogonal elastic 
modes are superimposed 

• Linear elasticity represented via a set of assumed modes (eigenvectors) and 
frequencies (eigenvalues) 

Other dynamic modes (propellant, actuator, and nozzle dynamics) are coupled to 
the base body 

• Complete system dynamics can be compactly derived from the Lagrangian 

The present formulation assumes 2 rigid DoF (rotation, translation), n slosh DoF, 
and /cflex DoF. 

• The system dynamics are implemented in first order descriptor form as 

E x — A x+B it , v— C x 4- D it 

• Controller, actuator, etc. are added via block diagonalization 

A vehicle with a single engine or cluster of engines can be analyzed as SISO at the 
gimbal break 
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Modeling Assumptions (III) 


• The aeroelastic model uses a quasi-rigid line load distribution 

Line loads derived from wind tunnel and CFD tests about zero angle of attack 

Does not account for nonstationary effects; e.g., the force distribution derivative is 
not a function of the flexible motion 

The continuous distribution is reduced to a discrete distribution integrated over 
intervals corresponding to FEM gridpoints attached to the load structure 
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Reference Frame Geometry 



elastic vehicle axis 


CM 




rigid vehicle axis 


System divided into empty vehicle with 
stationary propellant (v) , sloshing 
propellant (s), and nozzle (n) 
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Linear Dynamics Model 


Derivation of the coupled dynamics yields a block linear model of the form 
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For a rigid rocket, the angle of attack is invariant along the centerline and is 
given by 
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Local Angle of Attack 


• For the flexible vehicle, the local angle of attack along the vehicle centerline is 
a function of rotation, lateral velocity, and local gust velocity 
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Model Dynamics with 
Distributed Angie of Attack 


• Model dynamics for rigid rotation, translation, and bending without other 
generalized forces and coupling (thrust, slosh, etc.) 
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• Linear model is modified with the following aeroelastic effects: 
Rigid body aerodynamic damping 
Rigid-flex coupling through aerodynamics 
Flex-flex coupling through aerodynamics (deorthogonalization) 
Wind gust inputs 
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Rigid Body Aerodynamic Damping 


• Rigid damping coefficients reduce to angular rate and translational terms, 
such that 
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• If the rigid line loads integrate to the 6-DoF force and moment linearization, 
the remaining coefficients are identical 
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Aeroelastic-Rigid Body Coupling 


Let a*b^\a x b x a 2 b^ ... 
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Rigid coupling from flexibility via aero reduces to 
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Flexibility 


Define weighted dyadics of the eigenvectors, such that 
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Flexibility cross coupling collapses to delta stiffness and damping, and 
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Finally, the flex response decoupled from rigid-body, actuator, and slosh dynamics 
is given by 
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Gust Response 


• Gust velocity is assumed small with respect to forward velocity 

• A special case of the distributed gust 

Assume wind velocities are equal along the centerline 

The same rigid-body and flex coupling terms can be used for wind gust inputs 


r 1 


■ 1 

O' 



q S k r 



V 

q S k t 

B r= 

qS . A 

■ ■ ■ y 

V 



L ■ 


1 , 


16 


Energy | Environment | National Security | Health | Critical Infrastructure 



From Science to Solutions 


© Science Applications International Corporation. All rights reserved. SAIC and the SAIC logo are registered trademarks of Science Applications International Corporation in the U.S. and/or other countries. 



Applications 


• Ares I 

Two stage crew launch vehicle 

• Shuttle-derived five-segment solid propellant 
booster 

• Saturn-derived upper stage powered by 
LOX/LH2 J-2X 

- 56,000 Ibm payload capacity 

- >3M Ibf thrust at liftoff 

FS and US two-axis TVC control 
supplemented by roll control thrusters 

~10 minute ascent to LEO 
Dynamically similar to classic LV designs 

• Ares l-X 

First stage test vehicle with simulated upper 
stage 

-120 second flight to Mach 4+ 

- Flew successfully in October 2009 

• Other derivatives and study vehicles 


Launch Abort System 


Orion Crew Exploration Vehicle 
*7 (Crew Module / Service Module) 


Encapsulated Service 
Module Panels 
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J-2X Upper Stage Engine 
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Spacecraft Adapter 
Roll Control System 
Interstage 

Separation Plane 

Frustum 



Sth Segment Simulator 
and Avionics 

4- Segment Solid 
Rocket Motor 

(Space Shuttle) 

Aft Skirt (Space Shuttle) 
with BTMs & BOMs 

Nozzle and thrust 
vector control 



Thrust < 40K lb.- force 
First Stage separation 
begins +123 secs 
at 130 k feet 


Apogee: 1 50 kft 


USS/CM/LAS 

uncontrolled 

descent 


Booster, Parachutes 
and 

Recovery 


Range: 128 nm 
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Summary of Results 


• Augmented equations of motion implemented in FRACTAL (Frequency 
Response Analysis and Comparison Tool Assuming Linearity) m-code tool 

• Optional aeroelasticity included in standard production runs 

First stage control design (gains/filters, sensor location studies, architecture trades) 
Monte Carlo stability analysis 

- Parameter optimization (slosh damping) 

- In-flight load relief performance and coupled response simulation 
Sensitivity analysis 

• General aeroelasticity trends as expected for a generic launch vehicle model 

- High margin for flutter and divergence 

- Increase in rigid pitch (yaw) damping 

Aeroelastic damping of pitch (yaw) lateral bending modes 

- Improved control stability 

Increased vehicle transient bending deflections with distributed loads 
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Pole Freq 


Aeroelastic Loci 


• Typical variation 1 .32x max q (NASA SP-8003) 
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gain 


Coupled Autopilot Stability Margin Analysis 


Pitch Axis (Rigid Aerodynamics) 

Pitch Axis (Distributed Aerodynamics) 



GM: gain margin 
PM: phase margin 


phase 
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Rate Rotational Damping 


Rigid Body Damping 


Rigid Aerodynamic Rate Damping [pitch] 
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Rate Translational Coupling 



Coupled Dynamic Response Analysis 


FRACTAL gust response, distributed aero 
Scale Factor 20x / 40x / 20x (angle / flex / slosh ) 



vehicle station 
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Summary and Discussion 


• A planar first-order linear aeroelastic model has been incorporated into the 
stability analysis dynamics engine previously presented 

Evaluation of aeroelastic coupling with other dynamics (autopilot, propellant slosh, 
nozzle and actuator) improves confidence in system stability 

Transient simulation allows for the evaluation of autopilot performance in the 
presence of external disturbances, e.g., gusts, steady-state wind, thrust vector 
bias, etc., in the presence of distributed loads 

Aeroelastic stability can be predicted in a single tool using standard metrics to 
determine susceptibility to flutter and divergence in the presence of closed-loop 

control 

The present approach is fast and computationally efficient 

• Various other approaches exist for determining the coupled response; future 
work will compare these techniques 

Reduced models derived from coupled CFD-FEM analysis 
- Time-domain nonlinear simulation 

Identification of structural response via flight test programmed test inputs (PTI) 
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